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Abstract 

■ The operational meaning of spacetime fluctuations is discussed. Classical spacetime ge- 

. ometry can be viewed as encoding the relations between the motions of test particles in 

^ \ the geometry. By analogy, quantum fluctuations of spacetime geometry can be interpreted 



(N 



in terms of the fluctuations of these motions. Thus one can give meaning to spacetime 
fluctuations in terms of observables which describe the Brownian motion of test particles. 
We will first discuss some electromagnetic analogies, where quantum fiuctuations of the 



QQ ' electromagnetic field induce Brownian motion of test particles. We next discuss several 

o ■ 

explicit examples of Brownian motion caused by a fluctuating gravitational field. These 
Q , the fluctuating geometry, and fluctuations in the expansion parameter given by a Langevin 



version of the Raychaudhuri equation. The fluctuations in this parameter lead to varia- 
tions in the luminosity of sources. Other phenomena which can be linked to spacetime 



^ ! fluctuations are spectral line broadening and angular blurring of distant sources. 
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I. INTRODUCTION 



It is well known that classical spacetime geometry can be mapped by the motion 
of classical test particles which move along geodesies in the given geometry. Here test 
particle will be understood to include both massive particles, which move on timelike 
geodesies, and light rays moving on null geodesies. The fundamental quantity needed 
to characterize a spacetime geometry is the Riemann tensor, which in turn can be 
characterized by the phenomenon of geodesic deviation. 

A basic feature of a theory which combines quantum theory and gravitation is 
expected to be fluctuations of the spacetime geometry. These fluctuations can arise 
either from the quantum nature of the gravitational field itself (active fluctuations), or 
from quantum fluctuations of the matter stress tensor 0, Q, |^ (passive fluctuations). 
In general, both types of fluctuations can be present simultaneously 0, lESQ- 

The key question which we wish to address is how to give an operational signifi- 
cance to fluctuations of the spacetime geometry. Test particles will no longer follow 
fixed classical geodesies, but will rather undergo Brownian motion around a mean 
geodesic. This Brownian motion can be described by the quadratic fluctuations of 
some quantity characterizing the variation from the mean geodesies. Some explicit 
examples will be treated below. 



II. ELECTROMAGNETIC ANALOGIES 

In this section, we will discuss some examples of fluctuating forces of electromag- 
netic origin which are useful as analogies to the effects of the fluctuations of gravity. 
Specifically, we will examine the Brownian motion of a charged test particle coupled 
to a fluctuating electromagnetic field, radiation pressure fluctuations of a mirror, and 
Casimir force fluctuations. 



A. Classical Brownian Motion 

Consider a nonrelativistic charged particle of mass m and charge q in the presence 
of an electric field E. If we ignore magnetic forces, the particle's velocity satisfies 



where in general x = x(t). Here we will assume that the particle does not move very 
far on the time scales of interest, so we can ignore the time dependence of x. 





If the particle starts at rest at time t 



0, then we can write 




(2) 
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FIG. 1: A correlation function typical of classical fluctuations. The correlation function 
C(r) is non-negative, and decreases monotonically. The characteristic width of C(r) is the 
correlation time, Tc- 

Now suppose that the electric field undergoes fluctuations, so the mean trajectory 
of the particle is obtained by averaging Eq. 

(v) = ^ r(E(x,t))dt. (3) 
m Jo 

The fluctuations around the mean trajectory in the i-direction are described by 

(At;5 = ^ f f [(^,(x,ti) E,{^M)) - (^.(x,ti)) (^.(x,t2))] dt, dh . (4) 
JO Jo 

Thus, the dispersion in the i-component of velocity is given as a double time integral 
of an electric field correlation function. Let C denote this correlation function, and 
assume that it is a function of the time difference r = |ti — ^2! only, so we can write 

C(r) = (E,(x, ti) Ei(x, h)) - (Si(x, ti)) (Ei(x, h)) . (5) 

In this case, we may write Eq. (ji} as 

{Av^)=2^^ (\t-r)C{r)dT. (6) 
Jo 

The simplest possibility is classical or thermal-like fluctuations. Here the correla- 
tion function C is non-negative, and decays monotonically on a timescale of Tc, the 
correlation time, as illustrated in Fig. ^ In this case, performing the integration in 
Eq. (jni will yield a result of the form 

{AvJ) = a 4 C(0) t , (7) 
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where a is a co nstant of order one. This is the famihar random walk behavior where 
Wrms = \/{^vf) cx y/t. Note that here we are ignoring any damping effects, which will 
eventually stop the growth of v^rns- 



B. Quantum Fluctuations 

Now we wish to consider a model in which quantum, as opposed to classical fluc- 
tuation, drive the Brownian motion of the test particle. It is not clear that there 
is any effect associated with fluctuations of the quantized electromagnetic field in 
empty flat spacetime. However, there is a nontrivial effect produced by changes in 
these fluctuations, as will occur in the presence of a boundary. The simplest such 
boundary is a perfectly reflecting plane, and the resulting effects on Brownian motion 
of a charged particle were recently discussed in Ref. 0]. Here we will summarize the 
results of that paper. Let the reflecting plane be at z = 0, so we need to consider 
motion of the test particle in both the longitudinal {z) direction, and in a transerse 
(x) direction. 

The appropriate correlation functions are the renormalized two-point functions, 
which are the differences in the two-point functions with and without the plate, 

C.ir,z) = {E^M E.(x,n) = , (8) 



and 



C.iT,z) = (i?.(x,t') i?.(x,t")) = {Eyi^,t') Eyi^,t")) = - , , . (9) 



4z2)E 



Note that here (E) = 0. These correlation functions are plotted in Fig. |21 They are 
singular at r = 2 z, corresponding to the roundtrip light travel time to the plate, and 
can be either positive or negative. 

Despite the singularity, which can be attributed to the use of perfectly reflecting 
boundary conditions and a plate with a sharp boundary, the integral in Eq. (jH)) can 
be performed by an integration by parts approach. The results, for t^z are 

and 

(Ai;2)^ 1+ 1. (11) 

Unlike the case of classical noise, here the mean squared velocities do not grow in 
time. This is required by energy conservation, as there is no source of energy from 
which the particles can acquire kinetic energy. The fact that (Af 2) ^ reflects the 
fact that some energy is required to set up the system at t = 0, that is, a transient 
effect. The mathematical reason that we do not find growing mean squared velocities 
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FIG. 2: The renormalized electric field correlation functions near a perfectly reflecting 
mirror are plotted as functions of the time separation r for fixed distance z from the 
mirror. The solid line gives the correlation function for the longitudinal components, C^, 
and the dotted line that for the transverse component, C^. The dimensionless quantities 
Cz and Cx z^ are plotted. Both functions are singular at r = 2 z. 

is that the correlation functions have both positive and negative regions which tend 
to cancel one another so that 



Unlike a thermal state, a quantum state such as we are considering here is highly 
correlated, with subtle correlations and anticorrelations. 

It may come as a surprize that that the integral of Cz vanishes, as this appears to 
be a positive function. However, the singularity at r = 2^; effectively contributes a 
negative contribution when the integral is defined by integration by parts. A simple 
example of this is the following: 



The apparently positive function 1/x'^, when defined as a distribution, has a negative 
part at a; = 0. We might be able to remove the singularities in and Cz by a more 
realistic model of the plate, such one which includes dispersion and surface roughness. 






(13) 



5 



However, if this were done, then Cz{t) would become a finite function with a negative 
region near t = 2z. Similar behavior was found for the mean squared electric field 
near a plate when the singularity was smeared out by position fiuctuations 

The mean squared displacements in both the transverse and longitudinal directions 
can be found by integrating the Langevin equation, Eq. twice and then squaring. 
In Ref. jH, it is shown that the mean squared position fiuctuations are, in the t^z 
limit. 



\n{t/2z) 



and 



{Az') 
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It is of particular interest to note that (Ax^) < 0. This can only be understood if we 
account for the quantum nature of the test particles, and amounts to a reduction is 
the usual quantum position uncertainty of the particle. The results summarized in 
this subsection can also be generalized to the case of two parallel plates [Tfll |. 



C. Fluctuations of Radiation Pressure and Casimir Forces 

The model discussed in the previous subsection can be viewed as an analog model 
for active metric fiuctuations in that the charged particle is coupled directly to the 
quantized electric field. Here we wish to discuss some situations which are analogous 
to the passive metric fiuctuations. These arise when quantum fiuctuations of the 
electromagnetic field stress tensor cause force fiuctuations on uncharged material 
bodies. 

One example of this is the fiuctuations in radiation pressure when light in a co- 
herent state is refiected by a mirror. This effect was first discussed by Caves 11, ^ 



using an approach based upon fluctuations in the numbers pf photons striking the 
mirror. Caves' results were later rederived and extended using an approach based 
upon the quantum stress tensor. This approach requires an integration of a state- 
dependent, but singular cross term, which must be performed by an integration by 
parts procedure. Consider a free mirror of mass m and area A on which a linearly po- 
larized beam of light in a coherent state with energy density p shines perpendicularly 
for a time t. The resulting uncertainty in the mirror's velocity is 

(A^2^ = 4^t. (16) 

This uncertainty is likely to be a signiflcant source of noise in future generations of 
laser interferometry detectors of gravity waves. When it is detected experimentally, 
it will constitute an observation of quantum stress tensor fluctuations. 

Another phenomenon due to electromagnetic stress tensor fluctuations are the 
fluctuations of Casimir forces, which have been discussed by several authors ij, [l^ 



ITgI ITtI I. An example is the fluctuation of the Casimir- Polder force on an atom in 
the presence of a reflecting plate discussed in Ref. ^TQI. Consider an atom of mass 
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m and static polarizability a at a distance z from the plate. The transverse velocity 
fluctuations are found to be 




47 n^a'^ 
768 TT^m^zS 



(17) 



and 




3787 h^a^ 



(18) 



3840 TT^m'^z^ ■ 



These results are analogous to Eqs. (fTUI) and (fTT|) for the charged particle case. Again 
we can find suppression of the usual quantum uncertainty. In both cases, there are 
subtle correlations in the quantum fluctuations, which is a typical feature of quantum 
stress tensor fluctuations lisl. 



III. FLUCTUATING SPACETIME 

Now we turn to a discussion of a few of the phenomena which can arise when 
spacetime geometry fluctuates. 



A. A Loophole in the Classical Singularity Theorems 

The singularity theorems proven by Penrose, Hawking and others show that for- 
mation of a singularity is required in classical general relativity provided that certain 
energy conditions are obeyed. It has long been recognized that quantum matter 
fields can violate these classical energy conditions, for example in quantum states 
with negative local energy densities. This allows the possibihty of singularity avoid- 
ance in semiclassical gravity, where a classical gravitational field is coupled to the 
expectation value of a quantum matter stress tensor. 

It is less well known that even very small fluctuations in spacetime also create a 
distinct loophole in the singularity theorems . The reason for this is that the proofs 
of the theorems all rely upon focusing arguments. However, a bundle of geodesies 
will not come to an exact focus when there are fluctuations around a background 
geometry. This does not mean that minute fluctuations will necessarily prevent the 
formation of a singularity, but the classical arguments for singularity formation need 
to be revised. 



B. Lightcone Fluctuations 

One particularly striking effect of spacetime geometry fluctuations is that there is 
no longer a fixed lightcone. Recall that the lightcone plays a central role in classical 
relativity theory, as the boundary between events which can be causally connected 
(timelike or null separations) and those which cannot (spacelike separations). Once 
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the spacetime geometry undergoes fluctuation, no matter how small, this strict dis- 
tinction can no longer be maintained. 

Here we will summarize an ar iproach developed in Refs. H HI HI, HJ. (For 



other related discussions, see Ref |2J,|25|,|26[.) Let us consider a nearly fiat spacetime 
with small metric fluctuations, which can be described by a correlation function 
{h^u{x)hpx{x')) , where we assume that {h^^{x)) = 0. The metric fluctuations could be 
due either to gravitons in a nonclassical state (active fluctuations), or to fluctuations 
of a quantum matter stress tensor (passive fluctuations). Consider a source located 
at r = ro and a detector at r = ri . Then the mean flight time of light rays between 
the source and the detector is the classical result, Ar = ri — tq. However, there will 
be a root-mean-squared variation around this mean value of At, given by 



At=^, (19) 



where 



(al) = ^(Ar)2 r dr T dr' nVn^n^ {h^Ax)h,x{x')) . (20) 

Here is the tangent vector to the mean null geodesic separating the source and 
detector, and ai is the first order shift in the invariant interval a separating the events 
of emission and detection. 

Here At represents a mean time delay or advance relative to the classical propa- 
gation time Ar on the background spacetime. Thus, if we send a sequence of pulses 
between the source and detector, some will take longer than the classical time, but 
some will arrive sooner. This raises some interesting issues concerning causality. Nor- 
mally, signals sent outside of the lightcone could be used to send information into 
the past. On a fiat background, the order of emission and detection of a pulse trav- 
eling on a spacelike path can be changed by a Lorentz transformation. In the case of 
spacetime fluctuations produced by a bath of gravitons or by quantum matter fields 
in a nonvacuum state, there do not seem to be any causal paradoxes. The source of 
the fluctuations defines a preferred frame of reference, and effectively breaks Lorentz 
invariance. The case of lightcone fluctuations in the vacuum on a Minkowski space- 
time is more subtle, and will probably require a more complete quantum theory of 
gravity to be understood. 



C. Fluctuations of the Expansion: the Langevin-Raychaudhuri Equation 

The tendency of a gravitational fleld to act like a lens and to focus a bundle 
of geodesies is described by the Raychaudhuri equation. Consider a congruence of 
either timelike or null geodesies with affine parameter A and tangent vector fleld k^^. 
The Raychaudhuri equation gives the rate of change of the expansion 6 along the 
congruence to be p7l | 

JO 

— = -R^^k^k'' -ae^- a,,a^^ + u.^uo'^'' . (21) 
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Here R^,y is the Ricci tensor, o"^'^ is the shear and uj'^'^ is the vorticity of the congruence. 
The constant a = 1/2 for null geodesies, and a = 1/3 for timelike geodesies. The 
expansion parameter 6 is the logarithmic derivative of the cross sectional area A of 
the bundle: 

so that 6 < describes the case of converging geodesies. For ordinary matter obeying 
classical energy conditions, the Ricci tensor term acts to decrease 6, corresponding 
to the tendency of gravity to focus light rays. The Raychaudhuri equation plays a 
crucial role in the proofs of the classical singularity theorems. 

We are interested in interpreting this equation as a Langevin equation in which 
the Ricci tensor fluctuates. A more detailed account of this work is given in Ref. @|. 
For the purposes of this paper, we will further assume that the shear and vorticity of 
the congruence vanishes, and that the expansion remains sufficiently small that the 
6*^ term can also be ignored, so we can write 

^ = -i?^,pr, (23) 

The variance of the expansion can be expressed as a double integral of the Ricci 
tensor correlation function as 

((A^)^) = {9') - {9)' = dX dX' C^,„^(A, A') fc^(A)r (A) fc"(A')fc^(A') , (24) 

where 

Cf,uap{x,x') = {R^u{x)Rap{x')) - {R^u{x)) {Rap{x')) ■ (25) 

Expansion fluctuations are sensitive only to the Ricci part of the curvature, and 
hence to the passive fluctuations induced by the stress tensor fluctuations. Thus when 
both active and passive fluctuations are present, the expansion fluctuations provide 
a signature of the passive part. 

As discussed in Ref. j2^, the Ricci tensor correlation function has both state- 
independent and state-dependent singularities when x and x' are null separated. 
These singularities can removed by replacing the double line integration in Eq. ()24|) 
by integrations over a four- dimensional worldtube, corresponding to the spacetime 
volume occupied by a flnite wavepacket. One then needs an integration by parts 
procedure to deflne the resulting four- dimensional integrals. 

Next we need a physical interpretation of the result for A9. One such inter- 
pretation is a description of luminosity fluctuations. The focusing of a bundle of 
rays, where gravity acts as a converging lens, increases the apparent luminosity of 
a source. Conversely, defocusing (a decrease in 9) decreases the apparent luminos- 
ity. Thus fluctuations in 9 produce variations in the observed luminosity of a source. 
This is a well-known phenomenon in astronomy. Density fluctuations in the Earth's 
atmosphere lead to the twinkling of stars ("scintillation" in astronomer's terminol- 
ogy). In principle, quantum fluctuations of spacetime geometry could also produce 
scintillation of distant sources. 
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In general, the fractional variation in luminosity of a source, AL/L, is given by a 
double integral of a 9-9 correlation function: 

dt' dt" [( 9{t') 9{t") ) - {9{t')) {9{t"))] . (26) 




JO 



However, in some cases such as classical noise, the fractional luminosity fluctuations 
are directly related to the value of A9, 

(^) ^sA9, (27) 

\ Li / rms 

where s is the flight distance. 

There is a distinct physical effect which may be estimated from a calculation of 
A9. This is the degree of angular blurring of a distant source due to passive metric 
fluctuations. In Ref. jl^ , an argument is given which relates the root-mean-squared 
variation in angular position. Ay?, of a source to the variance of 9, with the result 

Av^rms = 2 A^rms • (28) 

This result is just a heuristic estimate, and a more precise formula for Aip will be 
derived in the next subsection. For another discussion of the propagation of light in 
a fluctuating geometry, see Hu and Shiokawa 0. 

In the Minkowski vacuum state, the result obtained from Eq. ()24|) is formally 
infinite, unless we average over a finite bundle of geodesies. Let this averaging be 
described by two length scales, a and h. Here a is the characteristic time during 
which the bundle of photons is emitted, and b is the typical cross sectional dimension 
of the bundle, that is, its size in a spatial direction perpendicular to the direction of 
propagation. The result is found to be of the form 

((A»)^. = ^J?|^. (29) 

where ^p is the Planck length, and is expected to be of order one. This leads to 
the estimates 

Even if we were to take a and h to be of the order of the photon wavelength, the 
smallest values they could reasonably have, this effect is too small to observe. It 
is plausible that one should never observe the effects of stress tensor fluctuations 
in the Minkowski vacuum, but that they should be masked by the usual quantum 
uncertainty of the test particles. 

A nonvacuum state, such as a thermal state, is a different story. In the case of a 
thermal bath at temperature T and a bundle of rays which are localized on a scale 
small compared to 1/T, one finds 

128v/q^ 



A9r^s = ^fpVTr^, (31) 

TT 
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where s is the flight distance and cq ~ 0.3468. Here we find the y/s behavior charac- 
teristic of a random walk. This leads to the estimate 



Although the effects of thermal stress tensor fluctuations are typically small, they are 
in principle observable far from the Planck scale. 

The possibility of observable effects of stress tensor fluctuations due to compact 
extra dimensions was discussed in Ref. j2^. The basic idea is that the smaller the 
compact extra dimension, the more violent will be the quantum stress tensor fluc- 
tuations of the Casimir energy due to the compactiflcation, potentially leading to 
observable angular blurring or luminosity fluctuations. It is true that if we hold the 
Newton's constant in the higher dimensional space flxed, one flnds that 

^drms OC %^ , (33) 

where G4+„ is the Newton's constant in 4 + n dimensions, and Vc is the volume of the 
n-dimensional compact subspace. However, in Kaluza-Klein theories, G'4+„ is related 
to the Newton's constant in four dimensions by 

G,+^ = VcG^. (34) 

This causes the observable effects to become independent of Vc, and in fact the same 
as in the Minkowski vacuum in a four dimensional spacetime. Thus one cannot detect 
the existence of compact extra dimensions through angular blurring or luminosity 
fluctuations. 



D. Redshift Fluctuations and Angular Blurring 

As has been noted above, two of the physical effects of the fluctuations of space- 
time are the angular blurring of images and fluctuations in gravitational redshifts or 
blueshifts. In this subsection, we will provide a unifled treatment of both phenomena 
in terms of the Riemann tensor correlation function. 

First consider a source and a detector moving along geodesies in a classical gravi- 
tational fleld, as illustrated in Fig. El Let be the four- velocity of the detector and 
A;^ be the tangent vector to the null geodesies separating the events of emission and 
detection. The detected frequency of photons will be proportional to — fc^t^. The 
constant of proportionality depends upon the choice of normalization (affine param- 
eter) for k^. Let 1 and 2 label two successive null geodesies, and fcf and ki^ be the 
values of k^ at the endpoints of each geodesic. Further choose the normalization so 
that k^ = (1,1,0,0) in the frame of the detector. This amounts to choosing the 
affine parameter to coincide with the detector's proper time at this point. Then the 
fractional change in detected frequency between 1 and 2 is 

[^)=-{k'2-K)t,. (35) 
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FIG. 3: A pair of null rays, 1 and 2 are emitted by the source S and detected by the detector 
D. Ray 1 is emitted at point A and detected at point while ray 2 is emitted at point 
D and detected at C . The tangent vector to ray 1 sX B \s k^, and that to ray S at C is 
The four- velocity of the detector is t^. Both the source and detector are moving along 
geodesic paths. The shift in the null vector is obtained as an integral of the Riemann 
tensor over the shaded region enclosed by the null and timelike geodesies. 

We are interested in the case where the frequency shift Auo/uj arises from the 
effects of gravity, rather than from any change in the output of the source. This 
can be enforced by requiring that the values of at the starting points of each null 
geodesic be related by parallel transport along the worldline of the source. Because 1 
and 2 are geodesies with tangent A;^, the values of k^ at the starting and end points 
of each geodesic are related by parallel transport. Thus if we parallel transport /cf 
backwards along i, along the worldline of the source to 2, and then along 2 to the 
detector, the result will be ^2 • Recall that if we parallel transport a vector around a 
closed path, the change in the vector can be expressed as an integral of the Riemann 
tensor over the area enclosed by the path. If the detector is in a flat spacetime region, 
then there is no effect from the parallel transport along the detector's worldline, and 
Afc^ = k2 — ki is the change after transport around the closed path. Let be 
defined off of the detector's worldline as the four-velocity of a congruence of timelike 
geodesies which include the worldlines of both the source and the detector. Then we 
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can write 

AF = J R^o^^p k" f da , (36) 

where the integration is over the 2-surface bounded by the four geodesic segments, 
the shaded region in Fig. IHl 

Now suppose that the Riemann tensor is subject to fluctuations (active or passive 
or both) described by the correlation function 

C^auf3 pyasix, X ) = {Rfiayfiix) Rpyasix )) — {RfMai^/six)) {Rp^yasix )) . (37) 

Then we can express the variance of the fractional redshift fluctuations as 

= ^ (^^) = J dada' C^^^pp^^s{x,x')t^'k"t''k'^ f'k^rk^ (38) 

where the indices fiauP are at point x and the indices p'ya6 are at x'. 

We can also relate the degree of angular blurring to the Riemann tensor correlation 
function. Let be a unit spacelike vector in a direction orthogonal to the direction 
of propagation of the null rays; thus s^t^j^ = s^k^ = 0. Then 

APs^ = tan V? ^ , (39) 

where if is the change in angle between 1 and 2 in the plane defined by the pair of 
spacelike vectors and = k'^ — f^. Here we are assuming that \ip\ <^ 1. We can 
express in terms of an integral of the Riemann tensor as 

(p = AFs^ = J R^c^up k"" k^ da . (40) 

The variance of due to Riemann tensor fluctuations can be expressed as 

A(^2 = (^^2^ _ (^^y = J^^ c^^^^ ^^^^^^^ 3,/) s'^k'^fk'^ sPk^rk^ . (41) 

In cases where the Riemann tensor fluctuation are dominated by Ricci tensor fluc- 
tuation, the above result can be expressed in terms of the Ricci tensor correlation 
function. In general, we expect the heuristic estimate Eq. to agree in order of 
magnitude, but not in detail. 



IV. SUMMARY AND DISCUSSION 

In this paper, we have discussed the idea that the operational meaning of stochas- 
tic spacetime lies in quantities which describe the Brownian motion of test particles, 
either massless or massive. Some electromagnetic analogies, such as the motion of 
an electron or an atom in the modified vacuum fluctuations near a reflecting plate 
provide useful insights into the types of effects one might expect in quantum gravity. 
We examined several specific examples of Brownian motion in a fluctuating gravita- 
tional field. These include lightcone fluctuations, whereby the flight time of pulses 
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between a source and a detector can fluctuate around the classical flight time. One 
of the striking features of this effect is that the rigid distinction between timelike 
and spacelike intervals can no longer be maintained when gravity fluctuates. We also 
reviewed recent work which treats the Raychaudhuri equation as a Langevin equation 
to compute fluctuations in the expansion of a bundle of geodesies. These fluctuations 
translate into luminosity fluctuations of a source. Thus distant objects would appear 
to "twinkle" when seen through a gravitational fleld with a fluctuating Ricci tensor. 
Finally, we developed a formalism which gives a unifled treatment of two other effects, 
redshift fluctuations and angular blurring. Both of these effects can arise when the 
Riemann tensor in the region between a source and a detector undergoes fluctuations. 

In most situations in the present day universe, these effects are small, but they can 
in principle be large far from the Planck scale. Furthermore, they could be signiflcant 
in the early universe and other strong gravitational fleld situations. 
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